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Abstract 

We study the betweenness centrality of fractal and non-fractal scale-free network models as well 
as real networks. We show that the correlation between degree and betweenness centrality C of 
nodes is much weaker in fractal network models compared to non-fractal models. We also show that 
nodes of both fractal and non-fractal scale-free networks have power law betweenness centrality 
distribution P{C) ~ C'^ . We find that for non- fractal scale- free networks 5 = 2, and for fractal 
scale-free networks 5 = 2 — l/d^, where ds is the dimension of the fractal network. We support 
these results by explicit calculations on four real networks: pharmaceutical firms (A^ = 6776), yeast 
(N = 1458), WWW {N = 2526), and a sample of Internet network at AS level {N = 20566), where 
N is the number of nodes in the largest connected component of a network. We also study the 
crossover phenomenon from fractal to non-fractal networks upon adding random edges to a fractal 
network. We show that the crossover length i* , separating fractal and non- fractal regimes, scales 
with dimension dB of the network as p~^/'^B ^ where p is the density of random edges added to the 
network. We find that the correlation between degree and betweenness centrality increases with p. 

PACS numbers: 89.75.Hc 
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I. INTRODUCTION 



Studies of complex networks have recently attracted much attention in diverse areas 
of science Many real-world complex systems can be usefully described in the 

language of networks or graphs, as sets of nodes connected by edges js, 6]. Although different 
in nature many networks are found to possess common properties. Many networks are known 



to have a "small-world" property 

flay. 

jetween any two nodes is very small. In addition, many real networks are scale-free (SF) 



10|: despite their large size, the shortest path 
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12j, well approximated by a power-law tail in degree distribution, P{k) ~ k~ 



where k is the number of links per node. 

Many networks, such as the WWW and biological networks have self-similar properties 
and are fractals 
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17| |. The box-counting algorithm [is], Q allows to calculate 



their fractal dimensions ds from the box-counting relation 



'B 



(1) 



where Nb is the minimum number of boxes of size Ib needed to cover the entire network 
(Appendix B). Structural analysis of fractal networks shows that the emergence of SF fractal 
networks is mainly due to disassortativity or repulsion between hubs 14|. That is, nodes of 
large degree (hubs) tend to connect to nodes of small degree, giving life to the paradigm "the 
rich get richer but at the expense of the poor." To incorporate this feature, a growth model 
of SF fractal networks that combines a renormalization growth approach with repulsion 



between hubs has been introduced [3]. It has also been noted 14] that the traditional 
measure of assortativity of networks, the Pearson coefficient r 191] does not distinguish 
between fractal and non-fractal network since it is not invariant under renormalization. 

Here, we study properties of fractal and non-fractal networks, including both models 
and real networks. We focus on one important characteristic of networks, the betweenness 
centrality (C), I20, 21, ; 



231 ] ■ defined as. 



(2) 



where (Jj^i) is the number of shortest paths between nodes j and k that pass node i and 
aj^k is the total number of shortest paths between nodes j and k. 

The betweenness centrality of a node is proportional to the number of shortest paths 
that go through it. Since transport is more efficient along shortest paths, nodes of high 



betweenness centrality C are important for transport. If they are blocked, transport becomes 
less efficient. On the other hand, if the capacitance of high C nodes is improved, transport 
becomes significantly better 24 1. 

Here we show that fractal networks possess much lower correlation between betweenness 
centrality and degree of a node compared to non-fractal networks. We find that in fractal 
networks even small degree nodes can have very large betweenness centrality while in non- 
fractal networks large betweenness centrality is mainly attributed to large degree nodes. We 
also show that the betweenness centrality distribution in SF fractal networks obeys a power 
law. We study the effect of adding random edges to fractal networks. We find that adding 
a small number of random edges to fractal networks significantly decreases the betweenness 
centrality of small degree nodes. However, adding random edges to non-fractal networks has 
a significantly smaller effect on the betweenness centrality. 

We also analyze the transition from fractal to non-fractal networks by adding random 
edges and show both analytically and numerically that there exists a crossover length i* 
such that for length scales i < i* the topology of the network is fractal while for i > i* it 
is non-fractal. The crossover length scales as £* ~ p~'^/<^B ^ where p is the number of random 
edges per node. We analyze seven SF model networks and four real networks. 

The four real networks we analyze are the network of pharmaceutical firms 25|] , an Inter- 



net network at the AS level obtained from the DIMES project 
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27|, PIN network of yeast 



30(1 . Pharmaceutical network 



28l . |29| and WWW network of University of Western Sydney 
is the network of nodes representing firms in the worldwide pharmaceutical industry and the 
links are collaborative agreements among them. The Internet network represents a sample 
of the internet structure at the Autonomous Systems(AS) level. The Protein Interaction 
Network (PIN) of yeast represents proteins as nodes and interactions between them as links 
between nodes. The WWW network of University of Western Sydney represents web pages 
(nodes) targeted by links from the uws.edu.au domain. Basic properties of the considered 
networks are summarized in Table. I. 

The manuscript is organized as follows: In section II, we study correlation between the 
betweenness centrality and degree of nodes, and we compare fractal and non-fractal networks. 
We analyze the betweenness centrality variance crc{k) of nodes of the same degree k and 
introduce a correlation coefficient R that describes the strength of betweenness centrality 
degree correlation. We also analyze the betweenness centrality distribution P{C) of several 
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model and real networks. In section III we study the transition from fractal to non-fractal 
networks with randomly added edges. Appendix A provides a short summary of the fractal 
growth model introduced in pj|. Appendix B discusses the box covering method and its 
approximations. 



II. BETWEENNESS CENTRALITY OF FRACTAL AND NON-FRACTAL NET- 
WORKS 

It is generally accepted [sij that in many networks nodes having a larger degree also 
have a larger betweenness centrality. Indeed, the larger the degree of a node, the larger the 
chance that many of the shortest paths will pass through this node; the chance of many 
shortest paths passing a low degree node is presumably small. Here we show that this is 
not the case for fractal SF networks. As seen in Fig. [Il^a) small degree nodes in fractal SF 
networks have a broad range of betweenness centrality values. The betweenness centrality 
of many small degree nodes can be comparable to that of the largest hubs of the network. 
For non-fractal networks, on the other hand, degree and betweenness centrality of nodes are 
strongly correlated. 

To demonstrate the difference in the relation between degree and betweenness centrality in 
real networks we compare original networks with their random (uncorrelated) counterparts. 
We construct the random counterpart network by rewiring the edges of the original network, 
yet preserving the degrees of the nodes and enforcing its connectivity. As a result we obtain 
a random network with the same degree distribution which is always non-fractal regardless 
of the original network. As seen in Fig. [Tt^b),the betweenness centrality-degree correlation 
of a random network obtained by rewiring edges of the WWW network is much stronger 
compared to that of the original network. Ranges of betweenness centrality values for a 
given degree decrease significantly as we randomly rewire edges of a fractal SF network. 

The quantitative description of the betweenness centrality - degree correlation can be 
given by the analysis of the betweenness centrality variance ac{k) attributed to nodes of 
specific degree k. Larger values of the variance (Tc{k) mean weaker correlations between 
degree and betweenness centrality of a node since nodes of the same degree have larger 
variations in betweenness centrality values. As seen in Fig. [21 in a region of small degree, 
betweenness centrality variance crc{k) of fractal networks is significantly bigger than that 
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of their respective randomly rewired counterparts which are not fractals. At the same time 
betweenness centrality variance of non-fractal networks is comparable or even smaller than 
that of the corresponding randomly rewired networks. Thus, the betweenness centrality 
of nodes of fractal networks is significantly less correlated with degree than in non-fractal 
networks. 

This can be understood as a result of the repulsion between hubs found in fractals fl^ : 
large degree nodes prefer to connect to nodes of small degree and not to each other. There- 
fore, the shortest path between two nodes must necessarily pass small degree nodes which 
are found at all scales of a network. Thus, in fractal networks small degree nodes have a 
broad range of values of betweenness centrality while in non-fractal networks nodes of small 
degree generally have small betweenness centrality. Betweenness centralities of small degree 
nodes in fractal networks significantly decrease after random rewiring since the rewired net- 
work is no longer fractal. On the other hand, centralities of nodes in non-fractal networks 
either do not change or increase after rewiring of edges. 

As seen in FigJU^b), the main difference in the betweenness centrality - degree correlation 
between fractal and non-fractal SF networks reveals itself in the dispersion of betweenness 
centrality values attributed to nodes of given degree, rather than in the average between- 



ness centrality values. 



36( 1 So, in order to characterize and quantify the overall betweenness 



centrality - degree correlation we propose a correlation dispersion coefficient: 

where (Jc{k) and cr^{k) are the betweenness centrality variances of the original and randomly 
rewired networks respectively and p{k) is the degree distribution of both networks. The 
dispersion coefficient R is the ratio between the mean variance < crc{k) > of the original 
network and < cr^(fc) >, the mean variance of the randomly rewired network. We note that 
fractal SF networks have bigger values of the betweenness centrality variance than their 
randomly rewired counterparts and therefore, have correlation dispersion coefficient R> 1. 
On the other hand (Jc{k) of the non-fractal SF networks is close or smaller than that of their 
random counterparts o"5(A;) which result in values of the correlation dispersion coefficient 
1 or i? < 1. The calculated values of the correlation dispersion coefficient R for the 
networks we considered in the paper are summarized in Table. I. 

The probability density function (pdf) of betweenness centrality has been studied for both 



Erdos Renyi [5|, l6|] and SF networks. It was found that for SF networks the betweenness 



centrahty distribution satisfies a power law 

P(C) ~ C-', (4) 



with typical values of 6 between 1 and 2 
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33, 



34( 1 ■ Our studies of the betweenness centrality 



distribution support these earher results (Fig. [3]). We find that 6 increases with dimension cIb 
of analyzed fractal networks. In the case of non-fractal networks, where ds = oo, estimated 
values of 6 seem to be close to 2. 

An analytic expression for 6 can be derived for SF fractal tree networks by using argu- 



ments similar to those used in 



34( 1 to find 6 for the minimum spanning tree (MST). Consider 



a fractal tree network of dimension ds- A small region of the network consisting of n nodes 



will have a typical diameter i{n) ~ n 



35| . Nodes in this region will be connected to 



the rest of the network via i{n) nodes. Thus, the betweenness centrality of those nodes is 
at least n. Since the number of regions of size n is N/n, the total number of nodes with 
betweenness centrality C > n in the network is 

4>(C > ~ i(n)— - n^^'^^-\ (5) 
n 

Thus, the number of links with betweenness centrality n is 

P(C) = A0~0(C + 1) -0(C) - C^/'^s-^. (6) 
Using Eq. (jll) we immediately obtain 

Thus, Eq. ([7]) shows that 6 increases with ds in agreement with Fig. [31 For non-fractal 
networks ds —>■ oo and 6 = 2. So non- fractal networks consist of relatively small number 
of central nodes and a large number of leaves connected to them. On the other hand in 
fractal networks, especially in those of small dimensionality, due to the repulsion between 
hubs, betweenness centrality is distributed among all nodes of a network. Analysis of the 
box covering method as a fractal test for some fractal and non-fractal networks studied here 
is shown in Fig. m 
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III. CROSSOVER SCALING IN FRACTAL NETWORKS 



We now study the behavior of fractal and non-fractal networks upon adding random 
edges. We analyze the crossover from fractal to non-fractal structure when random edges 
are added. To this end, we study the minimal number of boxes Nb of size is needed to cover 
the network as a function oi Ib as we add random edges to the network. Fig. [5](a) and [5]^b) 
show that the dimension oi the networks does not change. However, the network remains 
fractal with a power law regime Nb ~ ^"'^^ only at length scales d. below £*, a characteristic 
length which depends on p. For £ > £*, the network with added random edges behaves as 
non- fractal with exponential decay Nb ~ exp(— The crossover length separating 
the fractal and non-fractal regions decreases as we add more edges [see Figs. El^a) andO^b)]. 
We employ a scaling approach to deduce the functional dependence of the crossover length 
on the density of added shortcuts p. We propose for A''^ the scaling ansatz 



NB{d.p)-t{pr'-F{^), (8) 



where 



, M^^s M < 1 

F{u)^{ (9) 
exp(— m) m 1. 

With appropriate rescahng we can collapse all the values of NB{i,p) onto a single curve [see 
Figs. \5lc) and[S](d)]. The crossover length £*{p) exhibits a clear power law dependence on 
the density of random shortcuts [Fig. [S^e)], 

tip)r^p-\ (10) 

We next argue that asymptotically for large A^, 

r = l/dB. (11) 

When a fractal network with nodes and E edges has additional A <^ N random edges, the 
probability of a given node i to have a random link is = 2A/N. The mass of the cluster 
within a size £ in a fractal network is Mc ~ i"^^' ■ The probability of Mc{i) possessing a ran- 
dom edge is P = {2A/N)Mc- Thus, at distances I for which {2A /N)Mc ^ 1 we are in the 
fractal regime. On the other hand, large distances i for which {2A/N)Mc ^ 1 correspond 
to the non-fractal regime. Thus, the crossover length t corresponds to {2A/N)Mc{t) ~ 1, 



which imphes £* ~ ^-^/<^b ^ p-i/ds ^ where p = A/N. Note that the values measured 
for the two fractal networks, shown in Fig. EJ^e), r = 0.46 [ds = 1.9) and r = 0.39 {(Ib = 2.3), 
are slightly smaller then the expected asymptotic values, which we attribute as likely to be 
due to finite size effects. 

IV. DISCUSSION AND SUMMARY 

We have shown that node betweenness centrality and node degree are significantly less 
correlated in fractal SF networks compared to non-fractal SF networks due to the effect of 
repulsion between the hubs. Betweenness centrality distribution in SF networks obeys a 
power law P{C) ~ C^^. We derived an analytic expression for the betweenness centrality 
distribution exponent 6 = 2 — 1/dB for SF fractal trees. Hence, fractal networks with smaller 
dimension dB have more nodes with higher betweenness centrality compared to networks 
with larger ds- The transition from fractal to non- fractal behavior was studied by adding 
random edges to the fractal network. We observed a crossover from fractal to non-fractal 
regimes at a crossover length i*. We found both analytically and numerically that i* scales 
with density of random edges p as £* ~ with r = l/rf^. 
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APPENDIX A: A FRACTAL GROWTH MODEL 

I — I 

A growth model of fractal SF networks was first introduced by Song et al. In the 

nn 

core of the growth model lies the network renormalization technique |13l . |14| |: A network is 
covered with Nb boxes of size Ib- Subsequently, each of the boxes is replaced by a node to 
construct the renormalized network. The process is repeated until the network is reduced 
to a single node. The fractal growth model represents the inverse of this renormalization 
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process. The growth process is controlled by three parameters: n, m and e so that: 



N{t) = n*N{t-l) (Al) 
ki{t) = m*ki{t-l), (A2) 

where N{t) and fcj(t) are, respectively, the number of nodes of the network and degree of 
node i at time t. The parameter e is the probability of hub attraction e = E^^uhs/E. In the 
present study we limit our consideration to two distinct types of models: fractal (e = 0) 
and non-fractal (e = 1). At each growth step we run through all existing nodes. With 
probability e we increase the degree of a given node by attaching (m — l)ki{t — 1) new nodes 
(this corresponds to hub attraction). With probability 1 — e we grow (m — l)A;j(t — 1) — 1 nodes 
using remaining node to repel hubs. Thus, the entire growth process can be summarized as 
follows (see Fig. [6]): 

(1) Start with a single node 

(2) Connect (m — l)ki{t — 1) extra nodes to each node i to satisfy Eq. flA2p . With 
probability 1 — e use one of the new nodes to repel node i from the central node. 

(3) Attach the remaining number of nodes to the network randomly to satisfy Eq. flAip . 

(4) Repeat steps (2) and (3) for the desired number of generations g. 

The networks constructed in this way are SF with 

A=l + i2I!i. (A3) 

iogm 

Fractal networks have a finite dimension 

logn , , , 

For derivations of Eqs. fjA3l) and fjAil) see Ref . Q . 

Here we refer to network models using a set of numbers (g,n,m,e). For example, the set 
(4, 5, 3, 0) should read as a 4th generation {g = 4) fractal (e = 0) network with n = 5 and 
m = 3. According to the above growth process for this example, (4, 5,3,0), N = = 625, 
E = N - 1 = 624, A = 1 + logn/logm = 2.46, and c/b = logn/log2 = 2.32. 
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APPENDIX B: MODIFIED BOX COUNTING METHOD. 



The box counting method is used to calculate the minimum number of boxes Nb of size ^B 
needed to cover the entire network of nodes. The size of the box, f^, imposes a constraint 
on the number of nodes that can be covered: all nodes covered by the same box must be 
connected and the shortest path between any pair of nodes in the box should not exceed is- 
The most crucial and time-consuming part of the method is to find the minimum out of all 
possible combinations of boxes. In the present study we use an approximate method that 
allows to estimate the number of boxes rather fast. 

(1) Choose a random node (seed) on the network. 

(2) Mark a cluster of radius ^B centered on the chosen node. 

(3) Choose another seed on the unmarked part of the network. 

(4) Repeat steps (2) and (3) until the entire network is covered. The total number of seeds 
A^' is an estimate of the required number of boxes Nb- 

We stress that the estimated number of clusters A^' is always less than A'"^, the minimal 
number of boxes needed to cover the entire network. Indeed, the shortest path between any 
two seeds is greater then the size of the box Ib- Thus, a box cannot contain more than one 
seed, and in order to cover the whole network we need at least A^' boxes. 

Even though A^' is always less or equal to Nb-, the estimate may be good or poor based 
on the order we choose for the nodes. In order to improve the estimation we compute many 
times A^' (typically 100-1000) and choose the maximum of all A^'. 

Figures m^a) andlD^b) demonstrate the application of the modified box counting algorithm 
to several fractal and non- fractal networks. According to Eq. ([1]), dimensions of the fractal 
networks are obtained by calculating the slope of the Nb{(!-b) function in log-log format. The 
calculated dimensions are slightly underestimated due to a finite size effect of the analyzed 
networks. 

Figure HJ^c) represents (i^ as a function of the inverse number of generations g of the 
model. As number of generations g increases calculated dimension (Ib approaches the value 
given by Eq. (IA4p . 
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A similar algorithm was introduced in Ref . [16|] . The authors of this algorithm argue that 
it provides the same dimension of the network no matter how the boxes are chosen. In our 
algorithm we intend to estimate not only the dimension of the network but also the number 
of boxes. Thus, we are seeking the maximum A^' out of many realizations. 
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Network Name 


N 


E 


A 


ds 


Category 


Model lnf(7,4,2,l) " 


16384 


16383 


3.0 


N/A 


Non-Fractal 


Model 2nf(6,6,2,l) 


46656 


46655 


3.6 


N/A 


Non-Fractal 


Model 3nf(8,3,2,l) 


6561 


6560 


2.6 


N/A 


Non-Fractal 


Model lf(7,4,2,0) 


16384 


16383 


3.0 


2. 


Fractal 


Model 2f(6,6,2,0) 


46656 


46655 


3.6 


2.6 


Fractal 


Model 3f(8,3,2,0) 


6561 


6560 


2.6 


1.6 


Fractal 


SF Model 


2668 


3875 


2.5 


N/A 


Non-Fractal 


Uni West Sydney WWW 


2526 


4097 


2.2 


2.1 


Fractal 


Pharmaceutical [25] 


6776 


19801 


2.4 


N/A 


Non-Fractal 


Yeast [28] 


1458 


1948 


2.4 


4.2 


Fractal 


AS Internet [26] 


20556 


62920 


2.1 


N/A 


Non-Fractal 



"See Appendix A for abbreviation. 

TABLE I: Properties of the networks studied in the paper. Here is the number of nodes, 
E the number of edges, A the degree distribution exponent {P{k) ~ k~'^) and (Ib is the fractal 
dimension. The notation of model networks is explained in Appendix A. We consider only the 
largest connected cluster of the network if the original network is disconnected. 
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FIG. 1: (a) Betweenness centrality versus degree correlation profiles of fractal and non-fractal 
network models. Note the broader range of betweenness centrality values of small degree nodes 
of fractal network compared to that of the non-fractal network, (b) Betweenness centrality ver- 
sus degree correlation profiles of Uni Western Sydney WWW (fractal) network and its random 
counterpart. The randomly rewired network is non-fractal. Betweenness centrality and degree are 
correlated much stronger in nodes of the random rewired network. 
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FIG. 2: Betweenness centrality variance ac calculated for both original and rewired networks as 
a function of node degree k. Every point of the plot corresponds to the betweenness centrality 
variance calculated for nodes of the same degree k and normalized over the corresponding average 
betweenness centrality value < C > of the original network. Each of the plots includes the value 
of the betweenness centrality- degree correlation dispersion coefficient R, see Eq.([3|). Note, that 
small degree nodes of fractal networks: fractal model 3f(a), yeast(c) and Uni Western Sydney 
WWW(e) have significantly larger variance of betweenness centrality compared to their randomly 
rewired counterparts which are non-fractals. On the other hand, small degree nodes of the non- 
fractal networks: non-fractal model 3nf(b), pharmaceutical (d) and AS internet(f) have betweenness 
centrality variance comparable or even smaller than that of their randomly rewired counterparts. 
As a result -R > 1 for fractal networks and i? < 1 or i? ~ 1 for nonfractal networks. Thus, 
betweenness centrality - degree correlation is W^ker in fractal networks than in non-fractals. 



10° 




FIG. 3: Bctweenness centrality distributions of (a) Fractal model 3f and WWW network (frac- 
tal) and (b) non- fractal model 3nf and pharmaceutical network (non- fractal) . The data have 
been binned logarithmically. Both fractal and non-fractal networks exhibit a power-law range 
of betweenness centrality distribution consistent with P{C) ~ C"^. The measured betweenness 
centrality distributions (data points) are in good agreement with analytically obtained formula 
S = 2 — 1/dB represented by the straight lines. In non-fractal networks we expect 5-^2 since 
ds — > oo. 
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FIG. 4: Box-covering method applied to (a) models: If, 2f, Inf and 2nf and (b) real networks: 
WWW, Pharmaceutical, and AS Internet. The log-log plots of the number of boxes Nb needed 
to cover the network as a function of their size I b show clear "power-law" behavior for the fractal 
networks. The calculated dimensions are presented in Table I. (c) The calculated dimension of 
fractal model If for different generations g of the same fractal model network. Calculated value of 
dfi approaches the expected value (ds = 2) as the number of generations increases. 
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(a) Fractal model If 
d =1.9 
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(b) Fractal Model 2f 
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FIG. 5: (a,b) Box-covering analysis of fractal models If and 2f with added random edges. Networks 
remain fractal for length scales smaller than certain crossover length i*. Above i* the networks 
are no longer fractals. The crossover length i* becomes smaller as we add more edges. (c,d) Data 
collapse oi Nsi^iP) for the two fractal models. Appropriate rescaling Nb{^b) — ^ cl{p)Nb{^b/^{p)) 
allows to collapse all the values of Nb{1,p) onto a single curve, (e) The rescaling function b(p) = i* 
for fractal models If and 2f as a function of p shows a power law scaling of the crossover length 
i* ~ p^'^. Calculated exponents are ri = 0.46 and T2 = 0.39 respectively. Calculated values are 
slightly smaller than the expected values due to finite size effects. 
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FIG. 6: Constructing (a) non-fractal and (b) fractal networks with parameters n = 6, m = 3. 
The key difference between a fractal and non- fractal model is "repulsion between hubs" . In fractal 
networks nodes of large degree prefer to connect to nodes of small degree and not to nodes of large 
degree. 



19 



